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Abstract 

We  consider  a  class  of  dynamic  models  for  elastomers  involving  nonlinear  viscoelas¬ 
tic  (hysteresis)  as  well  as  nonlinear  finite  elastic  components  of  the  constitutive  laws. 
Existence  and  uniqueness  results  are  presented  along  with  sample  numerical  fits  to 
experimental  data  to  demonstrate  the  efficacy  of  the  models. 
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1  Introduction 


In  this  paper  we  present  a  well-posedness  result  for  a  new  model  describing  the  behavior 
of  nonlinear  elastomers.  These  models  for  filled  (e.g.,  with  carbon  black)  rubber- like  elas¬ 
tomers  are  motivated  by  formulations  for  large  deformation  tensile  dynamics.  When  properly 
modeled,  these  dynamics  require  finite  strains  [31,  33] 


du  1  /  du  \  2 

dx  2  \dx  J 


6  — 6 


(i.i) 


in  constitutive  laws,  where  e  =  0  is  the  usual  infinitesimal  strain  of  linear  elasticity  (see 
[8,  14]  for  basic  modeling  techniques  as  well  as  specific  examples  related  to  longitudinal 
(tensile)  deformation  models).  Computational  and  experimental  results  [9,  10,  11,  13]  have 
confirmed  that  nonlinear  constitutive  laws 


a{t)  =  Ge{e{t ))  (1.2) 

are  adequate  but  required  to  model  even  the  small  deformation  elastic  dynamics  for  lightly 
filled  elastomers.  However,  for  more  highly  filled  elastomers  (the  primary  focus  for  both 
passive  and  active  damping  devices),  hysteresis  is  present  and  viscoelastic  behavior  must  be 
combined  with  basic  nonlinear  elastic  behavior  in  constitutive  laws  of  the  form 

a(t)  =  Ge(e(t ))  +  Gv(e(t  +  s)  :  — oo  <  s  <  0).  (1.3) 

In  light  of  (1.1),  this  is  equivalent  to  attempting  to  determine  nonlinear  maps  gei  gv  so  that 
the  constitutive  law  is  given  by 

a(t)  =  ge(e(t))  +  gv(e(t  +  s)  :  -oo  <  s  <  0).  (1.4) 

We  choose  to  model  this  hysteresis  using  a  Boltzmann  type  law  of  the  form 

a(t)  =  ge(e(t))+  f  Y(t  -  s)-^-gv(e(s),  e(s))ds,  (1.5) 

where  a  denotes  the  stress  and  e  the  infinitesimal  strain.  This  stress-strain  law  implies 
that  the  stress  depends  not  only  on  the  current  strain  but  also  on  the  history  of  the  strain 
and  the  strain-rate.  This  law  contains  several  standard  internal  strain  or  internal  variable 
formulations  as  special  cases.  The  ADF  models  of  Lesieutre  [24,  25]  for  composite  materials 
exhibiting  both  elastic  and  anelastic  displacement  fields  are  formulated  on  the  assumption 
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that  the  host  elastic  material  contains  anelastic  materials  with  internal  strains  which  are 
elastic  strain  driven,  i.e.,  the  constitutive  laws  have  the  form 


a[t)  —  /'.i  (,(/•)  —  E2ei(t),  (1-6) 

where  the  internal  strain  is  described  by 

i^t) —  c2e(t) ,  (1-7) 


or 

ei(f)  =  f  c2e~Cl{t~s',e(s)ds. 

Jo 

More  generally,  one  might  have  a  nonlinear  version  of  (1.7)  given  by 

7|.(/-)  +  <V  i(0  = 

On  the  other  hand,  Johnson  et.  ah,  [22]  propose  a  linear  alternative  to  (1.7) 

f  i  ( /■ )  +  --  Cit(t) 


(1.8) 


(1.9) 


in  which  the  internal  strain  is  strain  rate  driven,  with  nonlinear  generalization  given  by 


£i(t)  +  c1e1(f)  —  gv(e(t)). 


(1.10) 


Our  attempts  to  use  models  equivalent  to  (1-6)  with  either  (1.8)  or  (1.10)  to  describe  either 
quasi-static  or  dynamic  tensile  test  data  [13]  have  not  been  successful.  However,  more  general 
models  of  the  form  (1.5)  which  correspond  to  an  internal  strain  model  of  the  form 

ei(f)  +  c1e1(f)  =  —gv(e(t),  e(t))  (1-11) 

have  performed  quite  well  in  modeling  our  quasi-static  and  dynamic  tensile  experimental 
data.  In  these  applications  we  required  that  the  derivative  ^  in  (1.5)  be  a  distributional 
derivative  (in  the  sense  described  below)  and  that  gv  have  the  form 


//«•(<(•“)*<(•“)) 


5„i(e(s)|  if  <(.s)  >0  1 
gvd(e(s))  if  e(s)  <0  J 


i.e.,  the  viscoelastic  response  differs  when  the  material  is  loading  from  that  when  it  is  unload¬ 
ing.  This  is  one  novel  feature  of  our  model  which  is  based  on  our  experimental  findings  [13]. 
In  the  following  we  shall  refer  to  the  nonlinearities  ge  and  gv  as  the  elastic  and  viscoelastic 
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response  functions,  respectively.  Our  kernel  function  Y  is  assumed  positive,  decreasing  and 
is  not  singular  at  0.  In  practice,  we  consider  y(£)  =  c2exp(  — $i^),  where  cl5c2  are  positive. 
Moreover,  we  expect  that  the  material  does  not  exhibit  infinite  memory,  but  depends  sig¬ 
nificantly  only  on  history  of  finite  length,  r.  Exponential  memory  kernels  Y  are  such  that 
y(£)  ~  0  for  /  >  r.  and  r  sufficiently  large.  Therefore,  we  approximate  (1.5)  by 

a(t)  =  //,(<(/.))  +  Y(t  -  s)-^-gv(e(s),  k(s))ds.  (1.12) 

In  a  typical  experiment  the  rubber  rod  loads  and  unloads  periodically.  Let  us  suppose  that 
it  loads  from  /, /<  to  tx+i  when  K  is  even,  where  t0  =  0. 

t 

i 

t  3 

t5 


O 

Figure  1:  Loading  and  unloading  of  the  rubber  rod 

Assuming  that  tx  <  t  <  tx+ i,  where  K  is  odd,  and  r  >  tmax,  where  tmax  is  the  maximum 
value  of  t  of  interest,  we  may  integrate  by  parts  in  the  integral  term  in  (1.12)  interval  by 
interval  to  obtain: 

v(t)  =  +  y(0)sw(eW)  +  /  Y{t  -  s)gv(e(s))ds 

Jt—r 

K 

+  E  Y(t  ~  h)(-l)k+1[gvi(e(tk))  -  gvd(e(tk))]-  (1.13) 

k= 1 

We  shall  refer  to  the  terms  in  the  summation  as  ’’jump  terms”  and  the  {//,.]  as  ’’turning 
points”.  We  note  that  we  obtain  a  similar  system  in  the  case  the  rubber  rod  is  loading,  i.e. , 
tic  <  t  <  l-K-jj  K  even. 
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We  consider  the  longitudinal  motion  of  an  elastomer  rod  of  length  £  and  let  x)  denote 
the  displacement  at  time  t  of  the  section  of  the  rod  originally  located  at  x,  0  <  x  <  £.  The 
motion  is  governed  by  the  equation  (for  a  careful  derivation  see  [8,  11]) 

pAcutt  =  Acax(x,t )  +  f(x,t),  (1-14) 


where  Ac  is  the  undeformed  cross  sectional  area,  and  p  is  the  density.  Combining  our  stress 
strain  law  with  (1-14)  we  obtain  the  model  for  /, /<  <  t  <  /. /x  _  | .  K  odd: 


pA, 


d2u  d 


dt2  dx 
K 


Ac(je 


du ' 
dx 


+  A, 


dx 


+AcJ2Y(t  -  tk)(-l)k+1[gvi  (|^(4) )  -gvd  (£(4)  )] 


Ar 


l  Bn  \  ft  l  Bn  \ 

Y(0)gvd  (  —  J  +  Ac^  ^  Y(t  -  s)gv  l  —  j  ds 

for  0  <  x  <  £  (1.15) 


9e  (toW)  + !  (0)9”i  (&w)  +  i-T(f  (to  w)  * 


u(t ,  0)  =  0 

n(0,  x)  =  9?0 

ut{ 0,  x)  =  0 

u(t,  x)  =  t  <  0. 


j  ds 

m 

(1.16) 

(1.17) 

(1.18) 

(1.19) 

(1.20) 

We  assume  that  the  rod  begins  its  motion  at  rest  with  possible  deformation  cy0 .  fixed  end  at 
x  =  0,  and  memory  <p1.  We  also  suppose  that  ge(C)  —  £  +  //,  (0:  311  (l  include  an  (internal) 
Kelvin- Voigt  damping  term  euxxt  (e  >  0)  in  the  model.  Then  the  system  in  variational  form 
is: 

PA°U  ~  £'4C““  “  (to  +  9'(  di*  +  5/(0)9*,(  to}  +  l-r  ^  “  SM  to(s))<b 


v/(/.,  0)  --  0  (1.22) 

u(0,  x)  =  (£>0  (1.23) 

ut(0,  x)  =  0  (1.24) 

u(t,x)  =  <p i,  t  <  0,  (1.25) 
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where  (1.21)  holds  in  the  sense  of  V *,  where  V  is  an  appropriately  chosen  Hilbert  space. 

In  the  following  sections  we  shall  show  that,  under  certain  assumptions,  the  system  (1.21)- 
(1.25)  has  a  unique  global  weak  solution.  The  study  of  Boltzmann  type  constitutive  laws  is, 
of  course,  not  new  and  similar  stress-strain  laws,  e.g., 

a(t.  :r)  -  -  (.’>(((/■■.  x))  •  f  a(t  —  s)tp(e{s,  x))ds  (1.26) 

J  —  OO 

along  with  the  resulting  equations  were  studied  by  numerous  other  authors.  To  summarize 
a  few  related  results,  we  note  that  it  is  known  that  in  case  the  kernel  a  is  regular,  globally 
defined  smooth  solutions  exist  for  sufficiently  smooth  and  small  data.  However,  smooth 
solutions  can  develop  singularities  in  finite  time  if  the  data  are  suitably  large  [15,  17,  18,  27, 
35].  The  existence  of  weak  solutions  was  shown  in  a  special  case  in  [29].  It  appears  that  if 
the  kernel  a  is  singular  at  zero,  more  regular  existence  properties  can  be  established.  (The 
use  of  singular  kernels  is  supported  by  several  investigators  in  polymer  physics,  e.g.,  Doi  and 
Edwards  [16],  Rouse  [32],  Zimm  [36],  and  Laun  [23].)  In  this  case,  existence  of  global  weak 
solutions  was  obtained  in  [2,  12,  19,  20,  28],  while  in  [21]  under  some  suitable  assumptions 
the  author  shows  that  global  smooth  solutions  exist  for  arbitrary  data.  A  good  survey  of 
results  concerning  viscoelastic  models  can  be  found  in  the  monograph  [34]. 


2  Existence  of  weak  solutions 


We  are  going  to  consider  the  system 

“  e“*‘  “  &  (&  +  +  5  <0)9,,(ar  +  i-A(<  “  S)9,WW)* 


+i;y(‘-4)(-i)i+iM^(4))-w(^(i4))])  =  m  ™  v-  (2.1) 

u(*,0)  =  0  (2.2) 

t/(a,  x )  =  (f0  (2.3) 

ut(a ,  x)  =  0  (2-4) 

u(t,  x)  =  < fjj  t  <  a,  (2-5) 


for  a  <  tx  <  t  <  tK+1  =  b.  Our  strategy  in  proving  the  existence  of  weak  solutions  is 
to  consider  the  intervals  [/.().,i|  ]•  ]/■  1 ,  /.  ^  ] .  ...,  [tK,ti<+ 1]  consecutively  (t0  =  0).  O11  [t0,B]  the 
rubber  is  loading,  gv  =  gvi  and  there  are  no  jump  terms  present.  First  we  are  going  to 
show  that  a  unique  weak  solution  1  *  exists  on  this  interval  and  is  smooth  enough  so  that 
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id 1  '(/-  I )..  - ( /.  | )  exist.  Next  we  suppose  that  a  unique  weak  solution  vdK\t)  exists  on  [0,  t^\ 

and  we  consider  the  interval  [ tx,tK+i ]■  There  we  have  a  similar  system  as  above  except  that 
now  we  pick  up  a  jump  term  and  p>0  and  p>i  are  modified,  i.e.,  (p0  =  :)  and 

f  y\  if  t  <  0  1 

^  I  u*A)  if  0  <  t  <  tji.  [ 


We  will  again  show  that  a  unique  weak  solution  u  exists  on  this  interval  with  the  necessary 
smoothness  and 


rdA)  if  t  <  t  k 
ii  if  I,  k  <  t  <  tf[+ 1- 


is  a  weak  solution  on  [0,  f/r+i]. 

We  choose  V  =  77^(0,  7)  =  {<f>  :  (j)  G  II 1  (0,  /),  </>(0)  =  0}  and  H  =  L2( 0,  t?) ,  so  we  have  the 
Gelfand  triple  V  H  V*.  (■,•)  denotes  the  inner  product  in  H ,  while  (-,-)v*,v  stands 
for  the  usual  duality  product.  Also,  there  exist  k.  k  such  that 


M\m  i2  <  n\Vv  <  ^ii^ii2 


(2.6) 


for  every  (j)  G  V.  Let  C[a,b\  denote  the  space  of  functions  u  :  [a,  b]  — >  H  such  that 

£[a,b\  =  {u-  U  e  L2(a,b]V)i  ut  G  L2(a,b]V)}. 

We  say  that  u  G  £[a,&]  is  a  weak  solution  of  (2.1)-(2.5)  with  (p0  G  V  and  ip1  G  L2[a  —  r,  a;  V) 
if  it  satisfies 


/  l  {">■.?>)  '  (  ^ X  )  T ‘-px)  T  {9ei^oc) i  T  (L  ( 0) gv (u^,  iix ) , 

J  a 

K 

+  (H  H(t  ~  tk)YiT  ~  tk)(-l)k+1[gvi{ux(tk))  -  gvd{ux(tk))\,  ipx)  dr 


k=i 


+  ( ut(t ),  <p>)  =  j  (F,  ip)v*}VdT  (2.7) 


+  (  /  Y(t  -  s)gv(ux,  ux)ds:  Lpx) 

Jr—r 

for  any  a  <  tx  <  t  <  tx+i  =  &  and  p>  G  jC [«..*]  as  well  as  the  initial  condition 

u(a,  x)  =  9?0- 


(2.8) 


Here 


H(s) 


1  if  s  >  0  1 

0  if  5  <  0.  J 
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2.1  Assumptions  and  Main  Result 


We  make  the  following  assumptions  (these  assumptions  on  the  nonlinear  term  are  the  same 
as  in  [6],  except  that  we  recpiire  that  a  strict  monotonicity  condition  hold  for  ge  +  gv.): 

Al)  ge ,  gvi  and  gV!i  are  continuous  nonlinear  mappings  of  real  gradient  type,  i.e.,  ge,  gvt .  gvfj  : 
H  — >  H  such  that  there  exist  continuous  Frechet  differentiable  nonlinear  functionals 


G’e,  Gvi,  Gvd  :  H  R  with 

G'Mf  =  R e(ge(<p),  ip)  (2.9) 

G’vi((f)4 ;  =  \lv(gr.,(^),t;:)  (2.10) 

G'vdW)i'  =  Re{gvd(<p)i  4>)  (2.H) 

for  every  ^  £  H. 

A2)  There  exist  constants  C ‘®,  C™,  C^d,  G|,  G”,  C%d,  G£,  G™,  C^d,  ve,  vm ,  vvd  >  0  such  that 

-fir1  -  *')IMI!  -  C\  <  G,(v=)  <  C2'M!  +  C*  (2.12) 

-  j(*_1  -  -  ef  <  <?„,(?)  <  cf M2  +  C3"  (2.13) 

-  j(*_1  -  OIMI2  -  cf  <  g„4v)  <  cflMI2  +  Cf  (2.14) 

A3)  The  nonlinear  functions  ge-,gvi,gvd  satisfy 

h'MW  <  CM  +  (2.16) 

IM^II  <  AlMI  +  A  (2.i6) 

IMv)ll<<51iMI  +  Ai  (2.17) 

for  every  £  H,  for  some  constants  G|j  C1^,  G^4,  G“,  (7^,  G^A 
A4)  The  Frechet  derivatives  of  ge ,  gv{  and  gv,i  exist  and  satisfy 

g'M  e  £(fr, ff)  with  <  <t,  (2.18) 

<&(¥>)  G  £(H,H)  with  ||^i(¥’)|U(H,H)  <  Gf  and  (2.19) 

e  with  \\gM\\c(H,H)  <  Cf.  (2.20) 

8 


A5)  We  suppose  that  the  following  monotonicity  conditions  hold  for  every  'tp  £  H ,  with 
\j  >  0: 

<*<'(//,  (Aj  -  geWiV  -  #)  +  4’)  >  /'Hr1  “  ^||2,  (2-21) 

(^)  -  //,  (</’),  <p  -  #)  +  Re(5w(<f)  -  //,,,:((/•):  F  «v)  >  /'Hr1  -  h’ll2-  (2.22) 

A6)  We  assume  that  V  is  a  smooth  memory  kernel  such  that  there  exist  constants  £74,  O5,  Ce 
such  that 

\Y(t)\<C4,  \Y(t)\<C5 ,  I  y(i)[<C6  (2.23) 

for  i  '  •  I  <-  tmax • 

A7)  The  forcing  term  A  £  T2(0,  tmax;  V*). 

Our  main  result  is  the  following: 

Theorem  2.1  Under  assumptions  Al)-A7)  there  exists  a  unique  global  weak  solution  of 
(2.1)-(2.5)  for  any  tp0  £  V,  £  L2(a  —  r,  a;  V). 

We  first  prove  that  under  these  assumptions  a  weak  solution  u  £  £[0jtl]  exists  for  any  <p0  £  V 
and  <pi  £  L2{—r ,  0;  V),  i.e. ,  we  use  a  —  t0  —  0,  b  —  tj  in  our  definition  of  the  weak  solution. 
The  method  of  the  proof  is  very  similar  to  the  one  in  [6].  First  we  give  an  apriori  estimate, 
then  introduce  Galerkin  approximations  to  the  weak  solution  and  justify  taking  the  limit. 
The  monotonicity  condition  (2.21)  plays  a  major  role  in  showing  that  the  limit  of  the  Galerkin 
approximates  is  a  weak  solution.  Next  we  show  that  this  solution  can  be  continued  on  [A,  f2], 
etc. 


3  The  weak  solution  on  [0,  t\] 

In  this  section  we  consider  the  interval  [0,ti],  where  the  rubber  is  loading  for  the  first  time, 
gv  =  fhi  and  no  jump  terms  are  present.  Thus  we  have  the  following  system: 

utt  ~  euxxt  -  —  +  ge(—)  +  V  (0 )<?«%)  +  V  (*  ~  s)gv(-^(s))ds 

=  F(t)  in  G* 
u(t7  0)  =  0 
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(3.24) 

(3.25) 


k(0,  x)  =  <p0 
Kt(0,  x)  —  0 
u(l.  x)  =  <pi,  t  <  0, 


(3.26) 

(3.27) 

(3.28) 


3.1  The  apriori  estimate 

Let  us  assume  that  1(0)  =  1.  We  choose  ip  =  ut  and  take  the  inner  product  of  (3.24)  with 

¥■ 


{ Utti  U t)v*,V  4“  (  :  Htx)  4“  Htx)  4“  (<7e(^a;)i  ^  '.%■/  4“  {Svi(Ux}  i  '^t:x ) 

4 if  1(4  s)gv{ux)ds,  —  (F[t),  ^ t)v*,v 


t—r 


Taking  the  real  part  and  using  (2. 9), (2. 10)  we  get: 


d 


dt  L2 


o  Ilwt(4)||2  4~  77 1 1  1 1 2  4-  Ge(ux)  4-  Gvi(ux) 


4“  ’U’tx) 


^ )  9v  (  5  ^tx)  ^ t)v*,v 


t—r 


(3.29) 


Integrating  from  0  to  L.  (where  t  <  ijJ,  and  using  (2.6),  (2. 12), (2. 13)  we  obtain 


IK(4)||2  +  v\\u\\l  +  2 e  f  ||ufe||2dr  +  2  /  Re(  /  F(t  -  <  M2 

Jo  Jo  Jo 

+M1||^o||y  +  2  f  Re(F,  ut)v*ydT  -2  f  Re(  f  Y(t  -  s)gvi(iplx)ds,  utx)dr ,  (3.30) 

JO  JO  Jr-r 

where  y  =  hL(ye  4-  vm).  Note  that  we  have  separated  the  hysteresis  integral  into  two  terms 
assuming  without  loss  of  generality  that  r  >  t.  Thus, 

||^t(4)||2  +  HMIv  4“  £  f  \\ut ||y^T  —  M2  4~  Mi  ||^o||2  —  2/j  +  2I2  —  2/g,  (3.31) 

Jo 


where 


hi  —  [  Rg(  f  1  (t  J ///  (  cr  1 .  uxt)dr, 

Jo  Jo 

/  R e(F,  ut)v*  vdr  and 
Jo 

f  Re(  f  1  (t  -^ ) *7 u ( ,  Utx)d 

JO  Jr—r 


h  = 


Ir  — 
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First  we  estimate  I1.  Let  M(t)  =  /0'  Y(t  —  s)gvi(ux)ds. 


\h\  = 


S  -^-vx)dr  <  (.w I  (y- M(r),ux)d 

Jo  dT  Jo  dr 


*  d 


+  ^||^a'||2  + 


(^(r)),  (/.,.(-))  +  {  /  Y(t  -  s)gvi(ux(s))ds ,  ux{t))  dr 


(3.32) 


Estimating  the  integral,  we  have 


-(Mr)),  (i.r( t))</t 


< 


||5«(^)IIII^IM7 


< 


|F(0)|  [\cr\\ux\\  +  Cf)\\ux\\dr  <  Ms  I"  \\ufydT  +  M4  (3.33) 

Jo  Jo 


and 


t  rr 


0  JO 


Y(t  -  s)gvi(ux(s))ds ,  ux(T))di 


t  rr 


< 


0  JO 


Y(t  -  s)gvi{ux)ds\\\\ux{T)\\di 


< 


/  f  C6(C?\\ux\\  +  Cll)ds\\ux\\dT  <  c6ct  ( /  IHI dr)  +  ever  /  II^IMr 

jo  jo  \Jo  /  jo 

<  M5  (  \\u\\ydT  +  M6.  (3.34) 

Jo 

We  use  a  standard  estimate  for  J2  : 

W  <  7  Jo  IM| vdT  +  ^JQ  \\F\\Ur, 

where  7  >  0  is  chosen  such  that  27  <  e.  Let  M(t)  =  /°_r  L’(r  —  •s)fi^((<£,if|)d.s.  We  have 

r=t  J  ,  8 


(3.35) 


\h\  = 


S  R ,e(M(T),utx)dT  <  ( M(t),ux )  j  (—M(T),ux)di 

Jo  T=0  Jo  at 


<  -^wmr + «kii2  +  ^ii^(°)ii2 + 4mm 


4  a 


4f 


+  /  (/  Y{t  -  s)gv((<f1)x)ds,ujj-  dr  <  f  ||  f  Y(t  -  s)gv((Lp1)x)ds\\ \\ux(T)\\d 

JQ  Jr—r  JQ  Jr—r 

+^mm2 + «iiwii2 + ^iiM(o)ir + mmm2  <  mtJ*  mut 


4/3 

+Ms  +  — 1|^(0||2  +  ak  1 \\u\\v  +  — ||M(0)||2 


(3.36) 
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Since 


is\\2  <  i  C5 


,  rcT\ 

Jo 


\u„ 


+  CJVr 


||M(t)||2=  ||  fY{t-s)gvi{ux)dU 

Jo 

<  Mg  f  || u || y dr  +  Mi0  and 
Jo 

\\M{t)\\  <Mn,  (3.37) 

by  (3.32)-(3.36),  we  have  that  (3.31)  yields 

|K(t)||2  +  i>|M|^  +  0  /  ||nt||ydr  <  C((^>0)  ^ivs^G)  +  f  ||n(r)||ydT  (3.38) 

jo  Jo 

for  a,  6,  7  small  enough,  (0  =  £  —  2y).  By  Gronwall’s  inequality,  it  follows  that 

|| (/,(/) ||2  •  h||»(/-)Hr  '  0  I  II "<|li  dr  <  C(rO:^i;  /fUi)-  (3.39) 

Jo 

3.2  Galerkin  approximations  and  their  convergence 

Let  ))f=1  C  V  be  an  orthonormal  system  in  V.  We  define  Galerkin  approximations  for 
(3.24)  by 

N 

uN{t)  =  J2ck  (*)$*, 

k=  1 


where  {c^(t)}  are  chosen  so  that  uN(t )  satisfies 


-^<UA  (*),  V’j)  +  {«*  ,  fax)  +  (&(«*)*  fax) 

■  (//r,('G  ):  VbO  +  {  /  -  s)gv(u^)ds,  =  (/•.  fa)y^ 

Jt—r 

for  j  =  1, .  .  .  ,  TV,  with  initial  conditions 

JV/n\  _  jV  ^  iV/n\  _  n 
ck  (0)  —  co*:  ~^ck  (0)  —  0- 


(3.40) 

(3.41) 


We  choose  {c^.}  so  that  lim/v_ot.  cm-'ll;k  =  Yo  in  V.  Arguing  in  the  usual  way  we  obtain 
that  the  Galerkin  approximations  also  satisfy  (3.39),  i.e., 


N 


(Oil2  •  '>|l"V(/)llv-  •  0  /  \\ut\\vd,T<C, 

J  0 


(3.42) 


where  C  is  independent  of  N . 
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It  follows  from  (3.42)  that  {uN}  is  bounded  in  (7(0,  /.],:•  I  )  C  /.2(0,  /. t ;  V  ),  and  {nf}  is 

bounded  in  (7(0,  /. t.;  // )  and  in  /.2(0.  /. |i  l  ).  Thus,  there  exists  a  subsequence  (denoted  again 

by  uA )  such  that 

uN  — *  u  weakly  in  L2(0,t1]V)  (3.43) 

uf  — >  ut  weakly  in  T2(0,  (i;  V).  (3.44) 

As  in  [6]  we  can  show  that  there  exist  a  subsequence  uN  and  u  £  Ao,*i]  such  that 

a)  The  set  {wA}  is  equicontinuous  and  bounded  in  (7(0,  ii;  V)  and 

uN(t )  — +  u(t )  weakly  in  V  (3.45) 

uniformly  in  t  £  [0,  /.  t  J .  i.e. ,  uN  — ►  u  in  C'n  (0,  '/j:  V ): 

b)  The  set  is  bounded  in  (7(0,  /.  j;  //  ),  equicontinuous  in  (7n  (0,  /.  i :  // )  and 

uf  (t)  — »  ut(t)  weakly  in  H  (3.46) 

uniformly  in  t  £  [0,fj]; 

c) 

— >  rq  strongly  in  L2( 0,ii;h/);  (3.47) 

d)  There  exist  he,  £  L2( 0,ti;  H)  such  that 

ge(u™)  — >  he  weakly  in  L2(01t1-1H)1  (3.48) 

flvii11*)  hvi  weakly  in  h2(0,  (•£  H).  (3.49) 

(3.48)  and  We  note  that  (3.49)  follow  from  (3.42),  (2.6)  and  (2.15),  (2.16).  We  can  prove 
(3.45)  exactly  as  in  [6]  using  a  general  version  of  the  Ascoli-Arzela  Theorem  [30].  The 
same  technique  can  be  used  to  show  the  statement  (3.46),  once  equicontinuity  of  {uA}  is 
established.  The  key  to  this  (as  in  [6])  is  the  boundedness  of  the  set  {uA}  in  T2(0,^;  V*)  : 
For  a  fixed  M .  let  =  YlkL l  ak{t)4’k^  ak  £  C^O,^].  For  N  >  M  we  have 

l“tt($Af)|  =  Jo  (utt(T)^M{r))V*,vdT  <  jo  l~(ux,  ®Mx)  ~  ®Mx) 

(//.(«*):  ®Mx)  -  (&»(«*),  $Mx)  -  { j  h'(r  -  s)9vi{Ux)ds,  $JV|S)  +  (/••  $m) 
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[0  1  ft  l 

-(y  Y(t  -  s)gv(iplx)ds,  dr  <  mn1  \\uN  (t)\\^\\^ M\\v dr 

JPt  1  /’tl 

!  ||itf  (r)||y||$M||y^T  +  rn3  /  ||^e(u^)|| ||$jw||yc?T  +  ra4  /  ||£f«(^)ll II$m|Mt 

0  Jo  Jo 

+m5  f  I  \\gvi{Ux)\\ds\\$M\\vdT  +  I  \\F(T)\\v*\\$M\\vdT  +  m6  f  \\<&M\\vdT 

Jo  Jo  Jo  Jo 

<  ii  H$Af  ||i2(0jtl.y)  (3.50) 

By  (3.42),  li  does  not  depend  on  N  or  M.  ((3.50)  is  valid  for  N  <  M  also,  by  the  orthogo¬ 
nality  of  'tpk-)  Equicontinuity  can  now  be  proved  exactly  as  in  [6]. 

Now  since  { u is  bounded  in  L2(0,ti;  V)  and  by  (3.50)  {u^t}  is  bounded  in  L2(0,G;  V*)., 
Aubin’s  lemma  ([6,  1,  26])  guarantees  that  the  set  is  relatively  compact  in  h2(0,t1;  H). 

Therefore  we  can  conclude  that  there  exists  a  subsequence  uN  such  that  uf1  — >  ut  strongly 
in  L2( 0,  ;  H),  which  proves  the  statement  (3.47).  Thus  the  statements  (3.43)-(3.49)  are  all 

established. 


3.3  The  limit  process 

Let  Vm  denote  the  class  of  functions  r/  G  £[0jtl],  which  can  be  represented  in  the  form 

M 

Tj{t)  =  ak{t)^k, 
k= 1 

where  a/.(t )  G  C1[0,ti].  Then  V  =  Um=1Vm  is  dense  in  £[0jtl],  By  definition  the  Galerkin 
approximations  u N  satisfy 

J0  |  («?  ■  n>)  +  (u",Tlx)  ■  £<«£:»/*)  +  (fh  («*):  V}  +  (<lvi{  "x  ): 

+{  f  Y(t  -  s)gm(u^)ds,rjx)  dr  +  (uf  (t),  g)  =  f  (F,g)v*ydT 
Jo  J  Jo 


Y(t  -  i!  //: 


(3.51) 


for  all  T]  G  Vm,  for  M  <  N.  Let  g  G  Vm  be  fixed,  with  M  <  N,  and  take  the  limit  in  (3.51) 
as  N  — »  oo.  We  obtain 

[  [~{ut,Vt)  +  {  ^xi  ^ ]x  )  T  Tjx^j  T  (/ie,  1]x)  T  ( hvi ,  Vx) 

Jo 

+  (  f  Y(t  -  s)hvids,Tfj  dr  •  (nt(t),g)  f  (F,  r])v*ydr 
Jo  J  JO 


Y(t  -  s)gv(Lplx)ds,gx)di 


(3.52) 
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We  can  argue  that  passage  to  the  limit  is  possible  in  the  hysteresis  term  since  for  given 

t  e  [0, t!],ipe  H ,  AT  :  L2( 0,  t-  //)  -►  tfh 


W(^)  = 


0  Jo 


y(r  —  s)4(s)ds,  ip)di 


is  a  bounded  linear  functional  on  L 2 f  0 .  /  :  //  )  and  therefore  weakly  continuous.  Consequently, 


E(r  -  r]x)d 


Y(t  -  s)hvi(s)ds,r]x)dT , 


JO  JO  JO  JO 

by  (3.49).  Passage  to  the  limit  in  all  the  other  terms  can  be  argued  based  on  the  convergences 
(3.43)-(3.49)  (for  details  see  [6]).  Equation  (3.52)  holds  for  all  rj  £  Vm  for  all  M,  and  since 
V  is  dense  in  £[o,*i],  for  all  rj  £  £[o,ti\  as  well. 

To  conclude  that  u  is  indeed  a  weak  solution,  it  remains  to  argue  that  for  any  rj  £  £[ 0jtl] 
and  for  t  £  [0,  tj] 


j  [(54  (  Ux)  7  Vx)  T  idvii  ux)i  hx) 

Jo 

+  (  /  Y(t  -  s)gvi(ux)ds7rix)  dr  =  f  [(he,  yx)  +  (h& rjx) 
Jo  J  Jo 


+  ( J 1  Y(t  -  s)hvids ,  //x) 
First,  we  prove  the  following  lemma. 


dr. 


(3.53) 


Lemma  3.1  There  exists  T0  with  tx  >  T0  >  0  such  that  for  any  t  <  T0  and  any  4,4  £ 
L2{0,h-,H)  we  have 


[  Re  [(//,  (d)  -  4  -  4)  +  ( c>)  -  <yr,(P),  4  -  4) 

Jo 


+  (  /  R'(r  -  a)(&40)  -  0  -  #)' 


dr  > 


f  jf 


dr. 


(3.54) 


T0  depends  only  on  ji  (from  the  monotonicity  assumption  (2.21)),  on  Cff  from  (2.19)  and 
the  bound  C5  on  Y . 

Proof:  By  (2.21)  we  have 

/  Re  )(//,  (d)  -  <1,  ftp),  4  -  4)  +  ( gvi(4 )  -  gvi{tf  ).  4  -  4) 

Jo 


+{ /  YiT  ~  s)(gvi{4)  ~  gVi{4))ds,  4~4) 


0 

t  j'T 


dr  > 


[  g\\4-4\\2dT 

Jo 


+ 


0  Jo 


Y(t  -  s)(gvi(4 )  -  gvi{4))ds,  4  -  4)di 


(3.55) 
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Now 


Y(t  -  s)(gvi(<f>(s))  -  gm(^(s)))ds,  0(t)  -  f  (r))dr  I 

/(/  Y(T~s)f  <lri{°(,>  '  0  ^)'/'j(^(*s)  V-(Y))(IO<Is,0(t)  f(-))</t 

«/0  */ 0  </0 


t  PT 


< 


o  JO 


|T(r  -  5)l  /  CTII^l5)  -  i>(s)\\d0ds\\(f>(T)  -  ^(r)||dr  <  C™C5  /  ||<^  -  V’lM' 


<  6™’Cst 


/  lk-^Hadr. 

Jo 


(3.56) 


Thus,  for  C^C5t  <  i.e.,  for  t  <  mWU,  -L  )  =  T0,  we  have  that  (3.54)  is  satisfied.  Thus 

°  ^  ‘65 ^  v  ^ 

Lemma  3.1  is  proved. 

Let  t  <  T0,  g  =  u N  in  (3.51)  and  take  the  limit  as  N  — >  00.  Using  the  convergences 
(3.43)-(3.49)  and  the  weak  lower  semicontinuity  of  norms  in  a  Hilbert  space,  the  comparison 
with  (3.52)  yields 


lim  Re  /  [<//,  (vr  ),  ur  ) 

A  — >oo  Jo  L 


+(dvi( u"),  u")  +  (l  Y(t  -  s)gm(u*)ds,  <v) 


N\i„  ,V\ 


dr 


< 


)  +  <//.  ;,  ux)  +  (  Y(t  -  s)hmds,  ux)di 


(3.67) 


Now  by  Lemma  3.1  with  <f>  —  7/'  and  if>  =  r/x  we  have  for  t  <  T0 

0  <  ^  Jo  ||«f  -  rjx\\ 2 dr 

~  ReJ0  “  9eM)i  ux  ~  lx) 

T  {{gvi{ux  )  —  (hi  (  l.r  ) ) :  ux  —  lx) 

+  (  /  T(r  -  s)^(u^ )  -  gvi{gx)ds,  u *  -  gx) 


dr. 


(3.58) 


Letting  N  — >  00  in  (3.58)  and  using  the  above  relation  (3.57)  along  with  the  standard 
Minty-Browder  technique,  we  find 


/  [(/A  (".r):  lx)  +  {9vt{ux),  gx) 
JO 


+  (  /  Y(t  -  s)gvi(ux)ds,  gx) 


dr 


dr  =  f  [(he,rjx)  +  (hvi,gx) 
Jo 


+  (  /  Y(t  -  s)hvids,riJ} 


(3.59) 
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for  all  rj  £  £[o,ti],  which  proves  (3.53)  for  t  £  [0,  T0],  We  note  that  it  is  possible  to  show  that 

in  h2(0,  T0;  if),  (3.60) 

by  setting  tj  =  u  in  relation  (3.58)  and  passing  to  the  limit  as  IV  — >  oo,  again  using  (3.57). 
Note  that  this  implies  ujv  — ^ ►  u  strongly  in  L2( 0,  T0;  V). 

Thus,  we  verified  the  existence  of  a  weak  solution  on  [0,  T0].  Since  T0  depends  only  on 
constants  from  our  assumptions  (and  independent  of  the  solution  itself)  we  can  extend  the 
weak  solution  to  [T0,2T0]  etc.  to  obtain  the  weak  solution  on  the  interval  [0,ti]. 


3.4  Uniqueness  of  the  weak  solution  on  [0,tj] 

Let  u  and  u  be  weak  solutions  of  (3.24)  on  [0,fi]  corresponding  to  the  data  <^0,  ^>l5 /.  Then 
w  —  u  —  u  satisfies  u;(0)  =  wt(0)  =  0,  w(t ,  x)  =  0  for  t  <  0  and 

("’«■.  V)v*,V  +  ("V:  '/.,••}  +  XI  Vx  )  +  (9e(Ux),  Vx) 

H-  {dvi  (^a;)  5  f]x)  ( J  Y (t  (^ar)  ds^  Tfx^j  ($e  3  T]x)  (dvi  (^a?)  5  Vx) 

~(J  Y(t  -  s)gvi(ux)ds,T]x)  =  0  (3.61) 

for  all  r]  £  jC[ 0jtl]  and  for  almost  all  t  £  [0,  tj],  by  the  dehnition  of  the  weak  solution.  Let 

A ge  =  ge(ux )  -  ge{ux) 

A gvi  —  gvi{ux)  gvi^y^x)  •  (3.62) 


d  ( 1 
dt 


Choosing  r/  =  wt  in  (3.61)  we  formally  (wt  £[0tl])  obtain: 

j^ll^tll2  +  ^ll^-||2|  +  t\\wtxf  +  (A  ge,  wtx)  +  (A  gvi,  wtx) 
+(Jo  Y(t  -  s)Agvids,  wtx)  =  0 
Integrating  from  0  to  L  ( t  <  1 1)  we  obtain: 

n«t(t)ir + *-ih«7(*)-||^+ a?  I'M* 

Jo 


(3.63) 


<  -2 


/[ 


(A ge,  wtx)  +  (A//,.,.  wtx)  +  ( j Y(t  -  s)Agvds,  wtx) 


dr 


(3.64) 
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Now  fg(Ag$j  wtx)dr  and  fg  (Agv ,  wtx) dr  can  be  estimated  as  in  [6],  i.e., 

1 


wtx )  dr 


o 

and 

t 

w>x)(h 

where  8  >  0  is  arbitrary.  Also, 

rt  rr  . 


< 


@S^~k  1  l  MrMdr  ■  ^  \\wtx\\2dT 

qfhi-1  J'  \\W(T)\\\,dT  +  c?s£  \\w,x\\2dT 


(3.65) 

(3.66) 

(3.67) 


Y(t  -  s)Agvids,  wtx)di 


to  JO 

rt  pr  ,  pi 


Y(t  -s)  g'vi{0ux(s)  +  (1  -  0)i7r(s))(w*  -  Ux)d0ds ,  wtx)di 
)  0  VO  Jo 

<  /  /  C,5^3’||u^|rfs||wte||rfT  <  •  h/-,  /  ||rt>te||2dr  (3.68) 

JO  JO  40  JO  JO 

for  any  8  >  0.  Hence,  if  we  choose  8  small  enough,  (3.64)  along  with  (3.65)-(3.68)  gives 

IKH2  +  &-1||u;(0||v  -  L  f  IMr)ll ydr.  (3.69) 

Jo 

By  Gron wall's  inequality  we  obtain  w(t)  =  0  on  [0,  t] ,  which  shows  that  the  weak  solution  is 
unique  on  [0,  ti]. 

Now  we  can  show  that  possesses  extra  smoothness,  i.e.,  u ^  G  (7(0,  /.j :  V”)  and  u j1*  G 
(7(0,  tj,  H).  We  consider  the  systems 


d 


w„  +  W  XX  +  ewxxt  =  F  -  —  (ge{u(xy)  +  y(0)fli„( 
+  ( Jt_r  Y(t  -  s)gv(u^\u^))ds^j 


iu(0,  x)  =  0 
u?t(0,  x)  —  0 


and 


~h  vxx  ~h  £-tvxxt  —  0 
w(0,  x)  =  c^i 
«i(0,  x)  =  0. 


(3.70) 

(3.71) 

(3.72) 

(3.73) 

(3.74) 

(3.75) 


18 


Since  £  (ge(u(V)  +  y(0)^(u£1),  uj.1')  +  St-r  Y(t  ~  s)gv(u£\  vS^)ds)  G  L2{ 0,H;  V *)  by  [7]  we 
get  that  a  unique  weak  solution  w  and  v  of  (3.  72)  and  (3.75),  respectively,  exists  on  [0,  ti] 
with  w.  v  G  (7(0,  If.  I-  )  and  wt,vt  G  (7(0,  H;!!).  However,  w  +  v  and  id1)  both  solve  (3.24)  , 
so  by  uniqueness  w  +  v  =  id1)  and  id1*  G  (7(0,  H;  V),  u G  (7(0,  H;  H ). 


4  Continuation  of  the  weak  solution 


We  suppose  that  a  unique  weak  solution  id K )  exists  on  [0,  /,/<]  with  id K ■  G  (7(0,  /./<:  V  ),  ujA  *  G 
(7(0,  / /,- :  17),  and  we  consider  the  next  time  interval,  [/, /<,  ] .  (We  note  that  the  jump  points 

tj  are  dehned  by  u(tj)  =  0,  hence  continuity  of  ut  across  these  points  is  guaranteed.)  We 
have  the  following  system: 


^tt  &'Ujxxt  +  Qei^x r)  ( 0 ) Qv  ( ^x 5  ^ x )  ~h  J  ^ 

+  y;(-i)i+iy((  -  4)M«i *>(<*))  - >(4))])  =  c(«)  in  v* 

jt=i  / 


i/(l,  0)  =  0 
i/(H,  x)  =  <do 
Kt(0,  x)  =  0 
i/((,  x)  —  (p  1,  /,  < 


(4.76) 

(4.77) 

(4.78) 

(4.79) 

(4.80) 


where  (^o  =  u(K\tn)  and 


<Pi 


<p  1  if  1  <  0 

idA)  if  0  <  t  <  I, k  ■ 


We  can  use  the  same  method  as  before  to  show  that  a  unique  weak  solution  exists  on 
[tic,  tic+i]  '■  first  we  give  an  apriori  estimate,  then  introduce  Galerkin  approximations  and 
show  that  they  converge  to  a  weak  solution.  If  we  examine  this  system  we  can  see  that  we 
can  perform  the  same  estimates  as  before  (except  that  now  we  integrate  from  /,  /<  to  /,  where 
tic  <  t  <  tx+ 1).  There  are  only  two  terms  that  are  different  from  those  in  the  previous 
calculations,  the  extra  jump  term  and  the  history  integral.  To  estimate  the  jump  term  we 
have: 


/'  ((-lG+WlT  -  >(<«-))  -  sMf- \tK))U,*)dT 

JtK 

<  /'  (-1)k+1F(t  -  tOIMnfhK))  -  «(»f)(‘A))llll«..IMn 

JtK 
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<CA\\gvi{u^\tK))-gvd{u{^\tK))\\l  \\utx\\dT<cj  \\utx\\2  dr  (4.81) 

JtK  JtK 

using  the  fact  that  riAK-  £  C(0,tK]  V )  and  the  assumption  A3).  The  history  integral  can  be 
split  into  three  terms: 

f  Y(t  -  s)gv(ux ,  iix)ds  =  [  Y(t  -  s)gv(tplx(s))ds  +  [  Y(t  -  s)gv(u^\ it(J^)ds 
Jt—r  Jt—r  JO 

+  f  Y(t  -  s)gv(ux,ux)ds.  (4.82) 

JtK 

Again,  all  these  terms  can  be  estimated  as  before  using  the  previous  bounds  on  tdAh  Thus 
we  obtain 

||wt(f)||2  +  +  $  f  \\ut Wv^t  <  C(ipo,  ipx,  (fen  (fj,  F,  f/c+ijf..  (4.83) 

JtK 

Galerkin  approximations  can  be  introduced  as  before  and  we  can  establish  the  same  kind  of 
convergences  as  (3.45)-(3.49)  now  in  the  interval  [//, .  /  /\  +  i].  The  analogue  of  Lemma  3.1  is 
straightforward  and  the  limit  process  can  be  carried  out  as  before.  Thus  we  obtain  a  weak 
solution  it  on  the  interval  [tK,tR+i],  which  is  again  unique.  We  can  argue  as  before  that  it 
has  the  required  smoothness  and  thus  the  proof  of  Theorem  2.1  is  complete. 


5  Numerical  Results 


The  model  mentioned  in  the  introduction  was  tested  on  a  series  of  physical  experiments  with 
filled  rubber  samples.  Here  we  report  on  one  type  of  dynamic  experiment  that  was  performed 
with  a  lightly  filled  rubber  sample.  Free  release  experiments  were  conducted  with  the  rubber 
rod  having  a  3  lb  tip  mass  at  one  end  (the  other  end  was  fixed).  Initially  the  sample  was 
lifted,  so  that  no  compression  or  extension  occurred,  then  the  support  was  removed  and 
the  mass  fell  freely.  (More  details  on  the  experiments  can  be  found  in  [10].)  To  model  this 
particular  sample  we  use  a  cubic  polynomial  for  <7e(e)  =  ape  +  a2e2  +  a3t3,  and  we  suppose 
that  gv  is  linear  and  does  not  depend  on  e,  i.e. ,  gv  =  gvt  —  gvfj.  Thus  we  use  the  following 
model: 


( f)u  \  ft  ( fin  \ 

+  AcY  (0)f/,;  (  —  j  +  Ac  V  ( t  —  s)gv  (  —  j  ds 


for  0  <  x  <  i 


(5.84) 


( 0u\  ft 

+  AfY  (Q)gv  (  +  4lc  J  V  (t  —  s)gv 


=  fj)  (5-85) 

X =i 
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u(t,  0)  =  0 

u(0.  x )  =  iq(0,  x)  =  0 
(/(/..  :r)  -  O',  t  <  0, 


(5.86) 

(6.87) 

(5.88) 


where  Y(t)  =  c2eClT .  Using  the  load  cell  data,  from  the  experiment,  we  set  up  a  parameter 
identification  problem  to  find  the  parameters:/?,  g1?  a2?  «s,  ci:  c2-  (Note  that  in  our  model 
(5.84)- (5.88)  we  did  not  include  Kelvin- Voigt  type  damping,  while  it  was  important  for  the 
theoretical  result.  We  can  think  of  its  coefficient  e  ~  0,  and  thus  reconcile  the  theoretical  and 
numerical  results.)  For  spatial  discretization  we  used  linear  splines,  while  piecewise  constant 
elements  were  used  in  the  time  discretization.  (More  details  on  the  computational  technique 
for  dealing  with  the  integral  term  can  be  found  in  [3,  4].)  We  used  Matlab  optimization 
routines  for  the  inverse  problem.  The  computed  result  shows  very  good  agreement  with  the 
collected  data,  (Figure  2). 


Figure  2:  3  lb  extra  weight,  relative  error  3% 

This  agreement  is  maintained  across  different  kinds  of  experiments  (quasi-static  as  well 
as  dynamic)  which  provides  evidence  that  this  type  of  model  works  well  in  approximating 
the  behavior  of  filled  rubber  samples.  A  paper  on  our  detailed  modeling  and  computational 
results  for  quasi-static  and  dynamic  experiments  is  currently  under  preparation. 
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